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Achievements 


OBiECTtVt-  1  -  HtGH-FlDELITY  CHARACTERIZATION 

•  Successful  Acquisition  of  MH/-P1V  Synchronized  with 
Near-  and  Far-Field  Acoustic  Data 

-  Ensembles  ol'  16  Single-Exposure  Images  with  t  fis  Time  Resolution 

•  Detailed  Jet  Characterization  through  Additional  Measurements 

-  Probe  Mapping  of  Pressure  Distribution  in  the  Plume 

-  Acoustic  Data  versus  Jet  0|xrraling  Condition 

-  Schlieren  How  Visualization 

-  Pulse-Burst  Laser  Flow  Visualization 


Demonstration  of 
Volumetric  P1V  with 
Plenoptic  Imaging 


Status: 

•  Dala  Reduction  of  Synchronous  Data  Set  in  Progress 

-  PIV  tillage  Processing 

-  Streak  Image  Processing 
-Acoustic  Data  Reduction 
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Overview  &  Objectives 
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Goal:  Enhance  the  timlerstanding  of  the  effect  that  near-nozzle  and 
inner-nozzle  tlow  conditions  have  on  jel  noise  radiation. 

OBJECTIVE  I  -  High-Fidelilv  Characterization  of  a  Heated, 

Over- Expanded  Supersonic  Jel. 

OBJECTtVL  2 -Source  Identification  Through  Developmenl  of 
Advanced  Analytical  Diagnostics. 

Objective  3  -  Enhanced  Computational  Modeling  of 
Hoi  Supersonic  Jets. 
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ACHIEVEMENTS 

objective  2  -  Source  Identification  Diagnostics 


•  Developed  “Shock  Deiection  Algorithm” 

-  For  application  to  the  pressure  signal  emitted  by  supersonic  jets. 

-  Enhanced  understanding  ot  ]x*enlial  sources  for  crackle. 

-  L n lut need  understanding  of  cumukilive  non-linear  distortion  -  s|xx*itic;illy  the 
implications  for  range  restricted environ ments  (lab  scale). 

•  Developed  Analysis  Methodologies  to  Investigate  Noise  Source  Terms 

-  Application  to  LES  data  shows  preliminary  results  that  wavenumtvr/lreqneney 
content  shows  consistency  along  lines  ol  constant  convection  velocity. 

-  Structured  with  the  goal  of  application  to  the  synchroniHis  data  set 
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ACHIEVEMENTS 


Objective  3  -  Enhanced  Computational  Modeling 

•  Enhanced  Hybrid  RANS/LES  Methodology  10  Include 
Inner-Nozzle  Features  in  the  Jet  Simulation 

•  Developed  and  Exercised  a  Computational  Framework  lor 
Phased  Array  Analysis 

-Also  applicable  iv  Objective  2. 
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Outline  of  Discussion 

•  Bnct  discussion  of  the  jet  facilities  and  characteristics 

-  Shock  Free  Nozzle  ai  UT-Austin 

-  Conic -Sec  I  ion  Nozzle  al  U.  Mississippi 

•  What  we  learned  about  Mach  wave  emission  and  crackle 

•  How  we  improved  Hybrid  RANS/LES  modeling  for  military  style  nozzles 

•  What  we  found  from  (he  computational  phased  array  analysis... 

•  What  we  have  learned  by  investigating  noise  source  terms  with  LES  data. 

•  How  we  setup  aud  acquired  the  synchronous  data  set. . . 

Preliminary  results  of  streak-image  analysis. 

-  iheliniimtry  rcsulis  from  Time* Resolved  P1V  analy  sis. 

•  Another  promising  approach  Plcnoptie  PIV 
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Shock  Free  Nozzle  at  UT-Austin 


Atm*  X  MOCwvzzif 

*  NPRl6  7l|.ilwck-fiw) 

■  ELsil  ditmuUa  bn 

*  Mass  flow  I  iMfyr  i 


tY.J^A.1. _ 

AhHUAt  llj  SU(  S 


V_  "  ■ 


2 


Shock  Free  Nozzle  at  UT-Austin 


Experimental  condition*  (e* hurried  from  quasi  l-D  isen tropic  compressible  flow  eq.) 
•  Current  data:  acquired  during  3  days  (Fall  2011) 

XrM-nrniy  (titty  1)  jykl*  array  ( titty  2) 
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Shock  Free  Nozzle  at  UT-Austin 

Spatial  topography  of  OASPL 

Allenlion  to  the  “fur-field  of  the  jet 

•  Classical  heart  shape  pattern  and  cone  of  silence 

•  Non-linear  metrics. 

•  Peak  noise  path  coincides  with  Mach  angle 

•  Skew  of  p  and  p 

•0  *0  SZv'T 

•  Kurtoiis  of  p  and  p' 

V  /  * 

•  Mortcy-Hou  ell  indicator 

•  Gol'dbcn  number 

J  «  (3 

•  Burgers’  equation 
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•  Cumulative  nonlinear  distortion  ' 
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Shock  Free  Nozzle  at  UT- Austin 


ifiE 


••  •'! 


:  OtiNKKr,...iv  KuJ^U^liKC.Vjui«Jkv>ivv*  Fn  -  12 


3 


2UI1  OSK  J<l  Noi«  HkC  UxtuJktmiW 

U.  Miss.  Anechoic  Jet  Laboratory 
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hr  r  irld  Arc  Arni.i  Mtcrupbunr* 

•  B&K  Type  4939  i/4-inch 
-  Free- Held  Microphones 

•  Type  2670  Pre- Amplifiers 

•  Curved  Rail  Support  for 
Ease  in  Positioning 
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Scar  Held  Microphone  Ftuillom 

»  Kulne  Sensors: 

-  Model  XT  140-  l(K)A  Transducers 

-  2.6  itim  diameter,  100  psi  range 
»  PCB  Sensors. 

-  Model  1 1 2A22  Accelerometer 
Compensated  Pressure  Probe 

-  50  psi  range,  Quartz  Sensing  Element 
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U.  MISS.  ANECHOIC  JET  LABORATORY 

SUMMARY  OF  TEST  CONDITIONS  CONSIDERED 

•  Jet  Nozzle 

-  Conic  Section  Converging-Diverging  Nozzle:  Design  Mach.  Mj  =  1.74. 

-  Exii  Diitineicr.  /)/  *  50.8  mm 

•  Pressure  Matched  Operation 

-  A  PR  =  5.2.  /TA*  =  3  37.  y  =  1 .372.  M,  =  1.74.  Ma  =  2.54 


•  Over-Expanded  Operaiion 

\RR  =  3  9,  TTR  =  3.37,  y  =  1 .368,  Mj  =  1 .55.  Ma  =  2.35 
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Mach  Wave  Emission  and  Crackle 

Challenges 

■  Crackling  and  non-crackling  signal  have  indistinguishable  spectral  content 
•  No  unique  measure  of  crackle  -  PDF  based 

•  Skewness  of  the  pressure,  Ffowcs  Williams  ct  at  1975  S(p)>  0.4,  S(p)  <0.3 

■  Skewness  of  the  pressure  derivative,  Mclnemy  19%,  Gee  et  al.  2007 

■  Kurtosis 


Objectives 

■  Isolate  and  study  the  shock  content  in  the  waveform 

■  Investigate  spatial  and  temporal  patterns  of  crackle 

■  Revisit  the  metrics  for  the  perception  of  crackle. 

■  *requires  a  robust  shock  detection  algorithm!  j 

n  i  .»  J  4  J 
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Shock  Quantification 


Wavelet  Transform 

■  Progressive,  complex-valued,  Morlet  wavelet: 

*(«//)  =  ,**/'«- H/«P/»  |*,|*6 

<'(M  =  (•>,?)”  l/-r  ('"■ 

■  Convolved  at  81  scales  (base-2  logarithmic  progression) 

•  Complex  valued  wavelet  coefficients: 


f'd 
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Energy  density  /  local  wavelet  power  speclrum  (WPS)  /  global  WPS: 

IC(I.I)  -f  E[f.i)  -S=l/TjrEi 
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Shock  Quantification 


Shock  Detection  Algorithm  (SDA) 


shock  rise  lime:  I"  f  1  mr.ioi 

shock  thickness,  A  III  /  i*.i*ln'  »< 


■  1.  pressure  gradient  \r  ;j.  m 

■  2.  user-defined  threshold  7  2.7o 

■  3.  median  time  of  p  and  p*  *  0 

■  Time-preserving  waveform  analysis  (wavelet  based)  to  study  spectral  properties 

{M.  «.«,  mm*  <4 
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_ _  Shock  Quantification _ 

Wavelet  Transform 

(■],  arc-array  mlc  5  (b) 
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Mach  Wave  Emission  and  Crackle 


Shock  properties  embedded  in  the  pressure  waveform  are  constant  along  rays 
emanating  from  the  post -potential  core  regions  of  the  How. 

Wavelet  based  transform  used  to  quantify  spectral  properties  of  the  shock  structure 

-  Increased  energy  in  higher  frequency  (nol  surprising  .  like  a  della  function1) 

Propose  the  use  of  percent  energy  gain  as  metric  for  the  perception  of  crackle 

-  Instead  of  skewness  ...  appears  shallow  angle  observer  bears  the  brunl  of  il  all. 

Wavepackel  models  are  resourceful  tools  that  demonstrate  a  real  promise  in 
understanding  these  mechanisms:  Morris  (2009),  Papamoschou  (201 1)  Kuo  et  al. 
(2012 ),  Jordan  &  Colonius  (2013) 


radiating  pressure 


pressure  standard  deviation 
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Quantifying  Crackle  -  Spectral  Character 


Scaling  to  replicate  full-scale 

•  Ffowcs  Williams  et  «/  1975.  "crackle  can  be  scaled" 

•  Scaling. 

n,  -  Im  C, 


Sl  m  llh. 

n'l>>  ~  (  o  C»*I n:Vt  tlMbu/o 
W.  *1  XT,*  1‘aaiA) 
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Energy  gain  -  krothapalli  et  al.  21)00 
GWSjrA  “>  Shock  spcclnun 
Only  function  of  0 

Ffowcs  Williams  ft  til  IV 74:  crackle  is 

mdeiMtntk'M  oj  the  distance  travelled  by  the  sound 
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Outline  of  Discussion 
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How  we  improved  Hybrid  RANS/LES  modeling  lor  military  slyle  nozzles... 
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CFD  Simulations  using  hybrid  RANS-LES 


Mi-an  and  Fluctuating  Chakacti  kistics 
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•  Exit  Boundary  Layer 

-Appears  initially  laminar 


•  Spcciral  Characterislics 

-  GiHKl  agreement  out  to  grid 
resolution  limit  of  sr  10  Ml/. 
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CFD  Simulations  using  hrles 

METHODOLOGY 
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•  Flow  Computation 

-  as  7  million  grid  points 

-  Computational  Domain: 

“Fan  I  ace"  lo  30  Dj 

5th  order  spalial  upw  ind  differencing 
-4,h  order  temporal  Runge-Kutia 

-  Parallelized  using  MP1 

•  Acoustic  Computation 

~  200  kHz  daia  recorded  on  an  acoustic 
data  surface  (ADS) 

•  CFD  lime  siep  *  5e-8  seconds 

-  Ffitwcs- Williams  and  I  lawk  mgs 
(FAN -I  I)  equations  used  It)  compuie 
noise  at  desired  fartield  location 
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Outline  of  Discussion 


t  Whal  wc  found  from  the  computational  phased  array  analysis. , . 
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Phased  Array  Calculations 


>  ADS  -j-  F-W&H  — »  Arriiy  Mies 

-  Possibility  to  use  locations  not 
feasible  in  a  laboratory  setting  due  to 
spatial  or  flowtield  restrictions. 

-  Allows  for  use  of  many  more 
microphones  than  practically  feasible 
in  experimental  work. 

-  Allows  for  rapidly  reconligurable 
phased  arrays  of  different  designs. 

•  Can  tailor  powtion/si/c  lo 

frequency/wavenunibcr/propagalion 
path  of  interest. 

>  Caveat; 

-  Obtaining  ’‘sufficient”  ensembles 
requires  numerous  time-resolved 
CFD  iterations. 


►  Application. 

-  Use  multi-arm  spiral  arrays  with 
sensor  at  logarithmically  separated 
spacing. 

Three  arrays  oriented  at  different 
polar  angles  to  localize  sources 
corresponding  to  different  jet  noise 
components. 
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Experimental  beamforming 

Estimate  of  Ensemble  Requiri  ments 


•  Data  Acquisilion: 

-  Linear  array  of  logarithmically 
separated  tensors  in  (.i\r)-plaiK 

•  Resulisof  Delay-aiiil-Sum  (DAS) 

Beam  forming 

-  DAS  results  relmed  using  DA  MAS 
technique. 

•  Conclusion 

-  Al  least  6  ensembles  at  A  f  =  -IS. S3 
Hz  are  required  for  satisfactory 
resolution  of  acoustic  sources 

•  Corresponds  lo  %2.5  million  CFD 
iterations. 


1 2S  rUMUilili  v.  \f  —  12.21 
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Experimental  bhamforming 

Estimate  of  Ensemble  Requirements 


•  Daia  Acquisition: 

Linear  array  of  logarithmically 
separated  sensors  in  plane 
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Insights  From  das  Beamforming 

Overexpanded  with  Centerbody  Hot  -  CFD 


•  General  dircciionulily  follow s  expectations. 

•  Centerbody  has  negligible  effect  al  these 
NPR  s. 

•  Shock  noise  clearly  present  in  the  overexpanded 
conditions  al  lhe  3.15  kHz  and  5  kl  1/  center 


frequencies. 

►  At  8  kl iz  all  three  configurations  have  similar 
noise  levels  across  all  frequencies 
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DAS  Results,  5  kHz  Center  Frequency 


•  Peak  noise  (135°):  mixing  noise  appears  to  peak  %  6 Dj 

•  Sideline  (90°) 

-  Peak  amplitude  «  8 Or 

-  Multiple  peaks  may  be  evidence  of  shock  cell  locations, 

•  Upsiream  (45°):  lower  amplitude  shock  noise  extends  up  to  «*  1 5Dj. 

©AiTROACWSTicS _ ^ _ HBS _ . 
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DAS  Results,  5  kHz  Center  Frequency 

Nozzle  with  Centerbody,  Pressure  Matched  Operation 


•  Peak  noise  (135°):  mixing  noise  appears  to  peak  as  6 D, 

•  Sideline  (90  ) 

-  Peak  amplitude  ss  8 Dj. 

-  Multiple  peaks  may  be  evidence  of  shock  cell  locations. 


•  Upsiream  (45  ):  lower  amplitude  shock  noise  extends  up  to  %  15 Dj 
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DAS  Results,  5  kHz  Center  Frequency 

Nozzle  with  Centerbody,  Pressure  Matched  Operation 


•  Peak  noise  (135°):  mixing  noise  appears  to  peak  «  6 Dj 

•  Sideline  (90  ) 

-  Peak  amplitude  «  8 Dj. 

-  Multiple  peaks  may  be  evidence  of  shock  cell  locations. 

•  Upstream  (45°):  lower  amplitude  shock  noise  extends  up  lo  Rs  15Dy. 

Q^SST’oTr^s  WM  4aa^.  . 
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DAS  Results,  3.15  kHz  Center  Frequency 

Nozzle  with  Centerbody,  Overexpanded  Operation 


•  Source  regions  extend  only  %  8 Dj. 

•  Shear  layer  &  shock  cells  seem  more  significant. 


-  Two  distinct  source  locations  in  the  shear  layers  apparent  at  4$ 


s-v 
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DAS  Results,  3. 1 5  kHz  Center  Frequency 


Nozzle  with  Centerbody,  Overexpanded  Operation 


•  Source  regions  extend  only  ~  8 Df. 

•  Shear  layer  &  shock  cells  seem  more  significant. 

-  Two  distinct  source  locations  in  the  shear  layers  apparent  at  45  , 


AfcftOA(.()lS]|lS 
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DAMAS  Results,  Nozzle  with  Centerbody 
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DAS  Results,  3,15  kHz  Center  Frequency 


Nozzle  with  Centerbody,  Overexpanded  Operation 


»  Source  regions  extend  only  %  8 Dj. 

»  Shear  layer  &  shock  cells  seem  more  significant. 

-  Two  distinct  source  locations  in  the  shear  layers  apparent  at  45  . 
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Outline  of  Discussion 


•  Whai  we  have  learned  by  investigating  noise  source  lerms  with  LES  data 


AmJALOUSULlL 


10 


itid  Rt.iuLlK®  RtU£*' 


Noise  Sources  in  Supersonic  Heated  Jets 


Challenges:  One  must  resolve  both  space  and  time. 

•  Experimental  handicaps: 

•  Difficult  to  resolve  both  space  and  time  (PIV-spatial,  LDV-temporal) 

•  L  imited  access  to  thermophysical  properties. 

•  Numerical  handicaps: 

•  Meshing  requirements  required  for  stabile  solution  (still  not  possible  for  DNS  at 
realistic  Reynolds  numbers). 

•  Statistical  convergence 

Objective  /  Motivation: 

•  How  we  can  leverage  the  strengths  of  both  tools  to  address  outstanding  challenges 
in  jet  noise9 

-  "L.ighthi 1 1-like” approach  (still  an  analogy!) 

•  Effects  of  viscosity  ,  heat,  compressibility,  2-D  vs  3-D 

-  Wavenumber  frequency  make-up  inside  the  source  Held. 


Vmi  (.mis  . 
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_  The  Large  Eddy  Simulation 

•  Fullv  expanded  tactual  aircraft  nozzle. 
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Large  Eddy  Simulation  provided  by  CRAFT-Tech 

•  Ideal  Gas  Law,  Sutherland  Law 

•  Time  Resolved  Data  Saved  at  200  kHz 

•  Differencing  scheme:  5th  order  in  space,  4lh  order  RK  in  time 

•  Spatial  deriv  ati\  es  2nd  order  -  central 


*uliO\K  At  K«hw  RcJuctk«i  bRC  Aiuuui  kcvww 

- Ejm 

Lighthill’s  acoustic  analogy 

•  1  nhomogeneous  wave  equation ; 

0-p 

■i  »  2p 

 i>2Tij 

in2 

°  ()Xi()xt 

1 

l 

v? 

The  Ligh thill  tensor:  = 


pu ,  it  j  pUyjj  -h  (p  \r,j\ 

(I)  (?)  (3) 


Viscous  stress  tensor  for  Stokesian  gas 


/  dui  duj  2  / duk\  ,  \ 

Ttj  *l  \  dxj  Oxi  3  dxk )  / 


How  do  each  of  these  terms  contribute  to  the  acoustic  analogy  and  the  generation  of 
sound  by  jets?  -given  the  restrictions  imposed  by  a  time-resolved  P1V  system 
(3)  viscous  effects? 

(2)  heating  effects? 


(I)  2D  versus  3D 


—  ill'll  j  TJ-’iU j 


■  “simple" 


prai  sbxzate*- 
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_ Radial  collapse  of  jet  properties _ 

•  Normalization:  r*  —  (r  -  nw)/ tow  —  no) 

•  Source  Convective  Speeds; 

-  Faster  than  the  mean  on  the  Low-Speed  Side 

-  Slower  than  the  mean  on  the  I  ligh-Speed  Side 

•  Horizontal  Lines;  Uv  =  a ^  Ur  =  0.7 Uj.  and  Ut.  =  0.8(7, 


FijiUrr  II.  (n)  ftmlitil  ill  trilxtl  l*ii  nf  Ilia'  niuvn'lUf  ti'Wil  V  nblmiH-al  from  Ihu-|xj1iiI  mrivlMl  Inim,  (l>)  C  ollupy* 
ut  t mcAii  vt'lorlly  itmlnlil  I  liv  nvurnui'  coiiv*tI  ivr  silmlly. 
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Wavenumber  Frequency  Characteristics 


Separating  radiating  and  evanescent  components  of  the  signal. 

-  Use  wavenumber  irequency  spectra. 

S(kxJ)  =  ^Z  I [  S(x,t)W{x)e-*k'x+2*Il)  dxdt 

-  Need  to  define  a  threshold  for  filtering1! 

•  Temperature  variations  effect  variations  in  the  sound  speed. 

•  What  sound  speed  do  we  use? 

-  That  of  the  ambient  field  or  that  of  the  izas  where  the  source  term  lies? 


<h*I  jit}  r  ** -  Th>  Mh'I  lint-  MnHllr>  (hr  k^nl  >*™l. 


A1KQ  \(  oi  Mil  s 


v. 


AlHMIJjSub. K^u.itlkL  \iuu.l km. 


Wavenumber/Frequency  -  Effect  of  Simplification 


ta  WxivwWmt  IhKimrr  nmtotir*  *r  (•)  •  (h>  <i  *  <r)  # i  (4)  Ai  *!«**  (ho  M<»  »W.  fV  *-«.  I«  * 
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_ Region  of  Acoustic  Efficiency _ 

•  Percentage  of  Acoustical  Efficiency:  ratio  of . . . 

-  Total  Speciral  Energy  (positive  ir  /): 

-  Acoustically  Efficient  Energy  sonic  threshold)  I  's 

•  Highest  Variance  in  Velocity:  Near  Center  of  Shear  Layer 

•  Highest  Efficiency  Percentage:  On  High-Speed  Side 


|Ih'  15  TW  i!kI<(  U  llw  |*t-ri<"iit<ii?€*  »«!>«■  i  wluk*  tU‘  l*  ft  <iiv*>  tin*  rwrgy  l**v*-l 

i»>tnuiliml  !*>•  It  nwixinitiKi  In  (  j)  mul  (<’  *t<*.  }/(l *,  —  m j.  )  in  |1>). 
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Present  Observations  from  CFD  Data  Analysis 


In  vvavenumber/frequency  space  it  appears  that  the  simplifications 
imposed  by  the  experiment  may  not  be  completely  limiting. 

LES  confirmed  negligible  effect  of  viscous  term  on  acoustic  analogy 
(Freund  2003). 

Present  analysis  firmly  sets  the  stage  for  investigating  the 
vvavenumber/frequency  spectra  using  the  time-resolved  PIV  data 
recently  acquired. 
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Setup  details  and  Run  Conditions 

1 

,  , 

•  Conic  without  Centerbody.  Overexpanded 

•  A7 >R  =  3.93.  Mj  =  1 .555.  Tj  =  1 350 °F.  Ma  =  2.34 

•  Cordiit  Acquisition: 

-303.  16- image  acquisition  1  every  15  sec.) 

-  1  fis  framing  rale  ( 1 5 /i>  total  lime) 

•  PCO.Edgc  Acquisition 

Image-Fair  (typical  IMV  method)  acquired  every  0.5  seconds 

-  First  Frame  -  Single  Exposure 

-  Second  f  rame  -  Multi-Exposure 

/JJS  XwmkJST imtwmM 

wAEMMt  <>fs||cs 

V  - 
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Streak-Image  Analysis 


Ra*«57 


IUM)- 


256-|)L\cI  ^iudoM  auiucoridatioii 


•  Application  of  window  ed  autocorrelation  methods. 

•  Identifies  direction  and  velocity  magnitude. 
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__ _ MHz  Rate  PIV _ 

Development  &  Implementation 


•  Utilize  a  MHz  rate  pulse  burst 
laser  system  and  MHz  rate 
camera  to  acquire  a  sequence  of 
particle  images  at  high  speeds 

•  Determine  velocity  using  PIV 
cross-correlation  algorithms 


•  Transport  equipment  form  AU 
to  NCPA 

•  Facility  integration  issues 

•  Dev  elopment  of  robust  image 
processing  algorithms  suited  for 
MHz  rate  image  data 
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_  _  MHz  Rate  PIV 

Challenge:  Particle  Seeding  &  Image  Quality 


P»*n01 


30 JJ  uK’k  Jtl  Noiw  koOuclKH.  Bkc  Annual  Ku\j«vi 

t'agw  faj 

MHz  Rate  PIV 

Dynamic  Evaluation  via  Ordinary  Least  Squares  (DEVOLS) 

•  MHz  rate  camera  image  quality 
compromised  by  particle  seed 
density,  laser  power,  intensifier 
and  beamsplitters 

•  A  Igorithm  de ve  loped  to  pro v  ide 
-  Robust  velocity  measurements 

llllllllllllllll 

-  Increased  dynamic  range 

•  DEVOLS 

*fc=i=J  T~ 

-  Utilizes  multiple  image  pairs 
centered  around  an  instant  in  time 
to  optimally  determine  local 
velocity 

-  Least  squares  fit  of  displacement 
vs.  time 

l  «unt  h|ohic i  lit  S\  nlhcliv  ' Noniy*  PIV  ikiU  uC  *n 

Itawn  vtNkn 

\  HLM 
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MHz  Rate  PIV 


Challenge:  Intense  Acoustic  Environment 

•  Pellicle  beam  splitters  used  to 
split  the  image  along  separate 
paths  have  a  resonant  frequency 
in  same  range  as  acoustic 
frequencies  produced  by  jet 

•  Led  to  image  blurring  due  to 
vibrating  beam  splitter 

•  Solution 

-  Construction  of  anechoic  chamber 
for  the  camera 

-  Move  camera  further  back  from 
flow  field 


ALM.W/QLMK'3 
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MHz  Rate  PIV 


Preliminary  Results:  DEVOLS  algorithm 
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MHz  Rate  PIV 


Preliminary  Results 


Streamwise  component  of 

velocity 

Standard  PIV 

Low-pass  filtered  in  space  and 
time 

Cropped  to  region  with  highest 
signal  levels 
Data  reduction  is  time- 
consuming  process  and  is 
ongoing 


K  2*  24  » 
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Outline  of  Discussion 


•  Another  promising  approach.  Plenoptic  PIV 


Kmimitimn  im  tkpm.4  3wMn 
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MHz  Rate  PIV 

Future  Work;  Better  Camera  Technology 

•  Cordin  -  222-4G 

•  Specialised  Imaging  -  Kirana 

-  Best  available  - 2  years  ago 

-  180  frames  (32  924  x  768  pixel 

-  16  frames  (0),  2,048  x  2  048  pixel 

resolution 

resolution 

-  Improved  image  quality 

—  Pellicle  Beamsplitters  -  frame-to- 

-  Custom  image  sensor  with  on-chip 

frame  jitter,  reduced  signal, 

memory  (no  beamsplitters  = 

susceptible  to  vibrations 

greater  sensitivity  and  no  frame-to- 

-  Intensifier  -  high-speed  gating 

frame  alignment  issues) 

decreased  spatial  resolution  and 

-  Available  for  demo  in  next  few 

image  quality 

months 

rilltl  *"  r*)Wc4Ac— Ac. 
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Volumetric  PIV 

Light-Field  Imaging  and  the  Plenoptic  Camera 

•  Analysis  of  planar  data  often 
must  assume  axisymmetric  or 
helical-type  flow  structures 

•  Volumetric  PIV 

-  Tomo-PIV- Multiple  cameras 

-  Generally  limited  to  thin  volumes 

!  (10  mm  or  less) 

•  Light-field  imaging 

-  Considers  the  complete  distribution 
of  light  in  space 

-  Record  both  the  position  and  angle 

‘  of  light  rays  traversing  a  volume 

•  Plenoptic  camera 

-  Utilizes  a  microlens  array  to  record 
the  light-field 

-  Small,  compact,  robust 

IV  ^ 

Plenoptic  Camera 

_ 'MM _ 
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_ Volumetric  PIV _ 

Computational  Refocusing  With  a  Flenoptic  Camera 
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SUMMARY 

•  Wluil  uc  Lamed  alxini  Mach  wave  cmivaoti  and  crackle 

•  Hove  we  improved  Hybrid  R  WS/LLS  nuxklmg  IW  nittit.tr>  tfyle  ik)//k‘>>. . 

•  NYlui  we  found  from  llx*  computational  plia\cd  array  uiutlysis. .. 

•  Wlial  we  ha\e  learned  hy  tnvcdiguiing  noise  xmv  lerms  with  LES  ilata... 

•  Strong  preliminary  results  from  Time  Resol  vet  t  PIN  analysis. . 

•  I  wiling  advancement  m  experimental  techniques. 


Sirwtti  Ctiln  hi  fkiywitl  At— tt\ 
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Volumetric  PIV 

Plenoptic  PIV:  Preliminary  Results 

•  Tomographic  reconstruction  + 

3D  cross-correlation  algorithm 

•  Proof-of-concept  at  NCPA 

-  Volume  size  60.7  mm  x  90  9  mm 
x  100  mm 

-  Centered  approximately  al  x/d  ■  J 

I  5 

•  1st  ever  volumetric  PIV 

measurement  in  a  supersonic  jet  1 
at  this  Reynolds  number 

•  Single  day  of  experiments 

-  Did  not  optimize  illumination, 
particle  seeding,  lomography 
parameters,  etc. 

-  Quite  pleased  with  results 
considering  this  was  I*  altempt 

O'  fm  B 

(  .1.11,  U  J. -  ■  1, .  ....  I/.JM.I 

(  ok.  conttpouk  lo  tilCJMOMtH  V.kKII) 
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Plans  for  Year  3 _ _ 

•  Reduce  synchronized  data  set  to  ensembles  of  time  dependent 
velocity  and  pressure  for  further  analysis. 

-  Apply  coupled  wavenumber  frequency  analysis 

-  Make  available  to  olher  BRC  teams  as  appropriale 

•  Modal  decomposition  of  velocity  data  retaining  short-time  evolution 
as  a  dependent  variable:  u,(x,\. r)  Li«<JV,(,)(v.y. r) 

•  Run  additional  CFD  iterations  to  further  enhance  computational 
database  for  beam-forming  analysis. 

-  Use  results  for  3-space  wavenumber/frequency  analysis. 

•  Perform  synchronized  data  acquisition  with  Kirana  camera 

-  No  later  lhan  Fall  2013  to  give  lime  for  analy  sis  prior  lo  end  of  program. 

•  Continue  analysis  of  Flenoptic  PIV  data:  Source  identification? 

•  Use  shock  detection  algorithm  with  stochastic  estimation  to  identify 
associated  flow- fie  Id  structures. 


iQ At  R()\(  ()1  M  kS _ BEJfl _ V. 


18 


Publications 


JU  (  "uNk'Ild  Raiuclkin  BRC  An 


*•*>  7J 


Jill  Baaks.  »  ) .  Tivvi  v,  C  LMiKKW.  S  »L,  Jaxms.  B  I  k  Pa\k  aa*.  p  <2oili  TW  ttil<xi  o*  Moa  oa  Tvrt>utr  «  Mua*  N»«w  w 
Sufv  rvjnW  ku,“  la  AIX\  \<nnp*n<  St  km,  r»  M<rtai*.  faivt  Mil  I0>  tOfUafc,  1 t- 

|2|  ItANHv  U  J  .TiWIV.C.  K  S  (2d|2l  V«ll»W4i  No.*  ProfM-alk*  ln«a  s  1-uHy  EvpumkJ  XL**  3  kf  Ui  UAkMAK 

t  Vtrmr  ittrUmg,  20|2  1177  (Saddle,  TN), 

in  MlkKU.N  LL1USS.C  W.TISNM  C  t  Dovaui.  B  Baaks  W  J  .Tmuo*  B  S.Mvvsts  It.  H  APaskuk  P.  (2UI2> 
“A  1  atodiiury  lur  Syiwtuuauu.  Skjii^m  htcU  Ak-ju  iin  m.1  Mil/  PtV  |r  Tcnfw/iMK  ttwtl  Cuiawn.  ktv"  * 

Prut  tttiutfc*  «ftht  tmrtnvtu-  201 2/ASSIt.  SXAO  !/«•(»;,  ASSIIAt  AD-IJ7II  Jim  iWl  1  N<x»  Y«*i  City.  SYi. 

HIKUIU.W  J  *Tl\MT  C.  fc  t20l2>~VJk»f  \l»i^  for  Noolinrj/ Supe«ijo»c  Vt  N«tw,“  ic  «airi«  «^/Ar  Amrrt.ai. 

*7  1 7.  Abarjtt  1)24  «  <S»i  Dwju,  CA  l 

|J|  IIMM1.K.  II  Bum  A.  B  A.  X  Till  KO».  B  S.  (2013)  "Ap$4K,  m»  at  VIII/  I  erne  Kj**  Hi  ■>  Dytun*.  R.  .-nje  IIV  te  >* 

Hijh  Trmfvntimr  Mm.  I  Conumnij  kt.  in  21"  ALKA  Am>»/v<  r  Sot*  ft  Mrrtmx.  Ml  34)774  (Orjpoux.  T\  i. 

|t»J  l*AXK  AH  P  LlOHvJ  SlMI  V.  N  ,  MlHK  VY  N.  I  ALMM1.C.U  <201 )»  Lo»  JtAatMrt  o<  A.ottOH  S.A*  *  M»*lk  t  oetm^tnt  In 
I'vmp  FNukJ  Ait/S  XtAMiarmralA,*  m  J/“  A  At  l  Asmyv.  r  Sntnt  fi  Afrrmt.  201  V**»l  3  iGr  irt  me.  T\v 

|7|  hair.  X  .Tivkiy.C  h.. Ml «h at. N.  t  Lvova. ow  a  foiu vk.  r <Mi 3i  Atuuai,  s<u.< lAkaon um ; u-a  « > i«U) 
fcvjjunJcd  «*1  Hc^U-d  Soivnomc  kt  u»  t*'  AilVCM-l  Arntanmu*  i  (  t^frnr;  <  r,  P/pri  201 3-2 193  i  Bc/^a,  Crntuny  i 

1*1  BaaH<  W  J  &TIVMY.C  C  <2013)  ^ujntdytai:  «r*Lfc  uidiiiinf  a.uu.1*:  tUViUMoccaardb)  *  fulty*ijuadni  NL*fc  3  jcC  io 

//  A/.UAI.U  Aw<tt.  i  c<«/r/rv<  201}  2l*M  (KcrUr .Orawttyl 

(9|  B.UKS  W  I  A  TlVNf  V.  C  E.<20l3t“A  T«npo«^a<idS^i»l(>uauiwMMMio4lWCn>k.t.|rCoaf>owat  MSupenoak.  JetNotw.'*  ia  2^ 

m  I'lliW.ifmMr'VW  ftu/niriir«i  W  Ownil.  p.  I/O 

|  |0J  Baakv  W  J  .20|V).A«,Hti(kt/rvKw//jjA  \vrJ  AtM  » ttk  I  ra,  IU.  Mil)  DiiViUtion. The  liinttMl)  d  Teta a  AiMm 

fin  l*AU.W.JnTwm.C.e  AKuKMN.M  S  -NO.IWJ.  rHOon**  of  AuvmMK  u  jsefonn*  ktat  ILfhSfKvA  irt*.'  JWwrf  .<  AW 

Sfftiwth  t,  |  in  nr*  K"»  | 


7/1/2013 


Qai  ko\<  ()l  Mil  s _ _ _ ^ ^ 


19 


